Global phase diagrams of frustrated quantum antiferromagnets in 
two dimensions: doubled Chern-Simons theory 



Cenke Xu and Subir Sachdev 
Department of Physics, Harvard University, Cambridge MA 02138, USA 

(Dated: January 2, 2009) 



Abstract 

We present a general approach to understanding the quantum phases and phase transitions of 
quantum antiferromagnets in two spatial dimensions. We begin with the simplest spin liquid 
state, the Z2 spin liquid, whose elementary excitations are spinons and visons, carrying Z2 electric 
and magnetic charges respectively. Their dynamics are expressed in terms of a doubled U(l) 
Chern Simons theory, which correctly captures the 'topological' order of the Z2 spin liquid state. 
We show that the same theory also yields a description of the variety of ordered phases obtained 
when one or more of the elementary excitations condense. Field theories for the transitions 
and multicritical points between these phases are obtained. We survey experimental results on 
antiferromagnets on the anisotropic triangular lattice, and make connections between their phase 
diagrams and our results. 
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FIG. 1: Distorted triangular lattice with nearest-neighbor Heisenberg exchanges J' (on all hori- 
zontal bonds) and J (on all other bonds), representing the geometry of systems examined in this 
paper. 

I. INTRODUCTION 

The study of exotic phases of quantum antiferromagnets has received a great impetus 
by the experimental discovery of a number of candidate S = 1/2 Mott insulators. The 
primary aim of our paper is to present an attempt to place the experimentally discovered 
phases in a single global phase diagram; such a phase diagram exposes new relations between 
the excitations of the various phases, and leads to theories for the possible quantum phase 
transitions between them. 

As will become clear from our analysis, we can generate distinct 'global' phase diagrams 
for distinct lattice types and exchange interactions in two spatial dimensions. We will present 
a general method for analyzing these, but will focus on a single lattice type, found in a num- 
ber of experimental systems: this is the distorted triangular lattice shown in Fig. dJ Thus 
we are interested in the S = 1/2 antiferromagnet, with SU(2) invariant Heisenberg exchange 
interactions J and J' illustrated in Fig. [TJ along with possibly additional longer range two- 
or multi-spin exchange interactions which have the same symmetry as Fig. [IJ A number 
of limiting cases of this lattice have been examined earlier, and we will connect with all of 
these results: 

(i) for J' ^ J the model becomes essentially equivalent to the square lattice antiferromag- 
nets considered in Refs. [iflsl, and our results agree with these earlier results in this regime; 

(ii) for J <C J' we have the quasi-one dimensional antiferromagnets which have been studied 
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in some detail by Starykh and Balents^ 5 -; 

(Hi) for J si J' we have the triangular lattice antiferromagnets for which our results will 
connect with those of Refs. 

(iv) Weihong, McKenzie and Singh— have performed a series expansion study for the entire 
range of J'/ J, and obtained phases which will also appear in our phase diagrams; 

(v) there have been a number of numerical studies^^ 2 -^ 1 ^ Q f isotropic triangular lattice 
case, J' = J, but with an additional four-spin ring exchange interaction, and our theory will 
provide candidate phase diagrams for this model; and 

(vi) the non- magnetic phases for J' = J have been modeled by the quantum dimer model 1 * 5 
on the triangular lattice by Moessner and Sondhi^, and our theory will also find their phases. 

Experimental examples extend over the full range of parameters for the lattice in Fig. 
and realize a variety of phases: 

• A remarkable series of experiments have been carried out by R. Kato and 
collaborators^i£»i2>2&2i'22^ on the organic Mott insulators X[Pd(dmit)2]2 (for a gen- 



eral review of the organic compounds, see Ref. |24j). Each site of the lattice in Fig. [T] 
has a pair of Pd(dmit)2 molecules carrying charge — e and spin S — 1/2. X ranges 
over a variety of monovalent cations, and the choice of different X allows experiments 
over a range of values of J' / J. The resulting phase diagram 2 ^ has magnetic order with 
decreasing critical temperatures from T c pa 42 K to T c pa 15 K across the compounds 
X = Me 4 P, Me 4 As, EtMe 3 As, Et 2 Me 2 P, Et 2 Me 2 As, and Me 4 Sb, as the value of J' /J 
increases from J'/ J pa 0.35 to J'/ J pa 0.7 (there are uncertainties in the overall scale 
of J'/ J, and there are also likely to be significant four-spin ring exchange terms). The 
magnetic order is likely of the two-sublattice Neel type^ 2 -, although there are no neu- 
tron scattering observations confirming this. The compound with X = EtMesSb has 
J' /J pa 0.85 has no observable Neel order— and has been suggested to be near the 
quantum critical point 22 at which the Neel order vanishes. Finally, the compound 2 ^ 2 ^ 
with X = EtMesP has J' /J ~ 1.05 has a ground state with a spin gap and sponta- 
neous columnar valence bond solid (VBS) order at low T. The VBS order vanishes at 
a phase transition observed to be at 25 K, and the low T spin gap is measured to be 
~ 40 K (the exchange constant J ~ 250 K). Thus these series of compounds appear to 
realize the Neel- VBS transition predicted in Ref. Q, and in this paper we will place this 
transition in the context of global phase diagrams of models on the lattice of Fig. [TJ 
We also note that a Neel- VBS transition as a function of increasing J'/ J has also been 
found in the series expansion study^. 

A separate set of experiments have been performed on the organic Mott insulators 
k-(ET) 2 Z, which also realize the S = 1/2 antiferromagnet on the lattice in Fig. [lj 
The compound with Z = Cu[N(CN) 2 ]Cl has J' / J pa 0.5, and has a ground state with 
Neel order—, as found in the Pd(dmit) 2 series above for small J' /J. The organic 
insulator with Z = Cu 2 (CN) 3 has J'/ J ~ 1, and appears to have no antiferromagnetic 



3 



or VBS ordering down to the lowest observed temperatures^ 1 ^ 1 ^ 1 ^. This is therefore 
a candidate for a spin liquid ground state, whose nature has been the subject of recent 
work— The bosonic Z 2 spin liquid state proposed for this compound in Ref. [32 



will appear in our phase diagrams below, and indeed will be natural point of departure 
for our entire analysis. We believe the experimental observation, noted above, of the 
other phases in our phase diagram can be regarded as a point of support for our 
perspective. We will briefly mention below how Z 2 spin liquid states with fermionic 
spinons^ -, and other related states, can appear in our approach. 

• The transition metal insulator—^ CS2CUCI4 has 5 = 1/2 Cu ions on the vertices of 
the triangular lattice in Fig. [I] with J m J'/ 3. The ground state has spiral antiferro- 
magnetic order, similar to that present in the perfect triangular lattice (J = J'), and 
as will appear in the phase diagrams below. An approach starting from the quasi one- 
dimensional limit J <^ J' has been successfully used^ to describe the spiral ground 
state, and also the inelastic neutron scattering spectrum at high energies. 

As noted above, our point of departure is a Z 2 spin liquid state. The earliest proposals 
of such liquids involved BCS-like states of paired, charge 0, S = 1/2 particles ('spinons') 
which were either bosons^ or fermions 3 -^. Fluctuations about this state are expressed in 
terms of a Z 2 gauge theory, in which the spinons carry a Z 2 electric charge, and hence the 
name of the spin liquid. A large number of other models of Z 2 spin liquids have appeared 
s i nce i£i^i^i^i39|4Sidi We will find it convenient to begin with Z 2 spin liquid in which the 
elementary spinons are bosons because it is connected naturally to a variety of ordered states 
found experimentally (which we have described above). We will denote the bosonic spinons 
by a complex field z a , where a =f, j is a spin index. 

Apart from the spinon, the other fundamental elementary excitation of a Z 2 spin liquid 
is a charge particle carrying Z 2 magnetic flux. This particle was pointed out in Ref. @, 
but its particular importance to the physical properties of Z 2 spin liquids was emphasized 
by Senthil and Fisher 36 , who called it a 'vison'. In all cases we shall consider, it is possible 
to combine the real visons into complex scalar fields v a , where a = 1 ... N v is an additional 
flavor index which depends upon the nature of the underlying lattice. The visons are bosons, 
but the spinons and visons have mutual semionic statistics 3 ^ 1 ^. Consequently, by forming 
bound states of the bosonic spinons, z a , and the visons, v a , we obtain S — 1/2 spinons which 
are fermions^,. This bound state formation offers a route to extending our analysis to the 
case of fermionic spinons, but we shall not comment further on this in the present paper. 

Our starting point is an effective field theory for the spinons, z a , and the visons, v a , 
which implements their mutual semionic statistics. As discussed generally by Freedman et 
alM, and implemented more specifically in Z 2 spin liquids with fermionic spinons by Kou, 
Levin, and Wen^ 3 -, this mutual statistics can be realized by a doubled U(l) Chern-Simons 
(CS) theory. A similar formalism was also applied to the cuprates, with a mutual statistics 
between spin and charge degrees of freedom^. To this end, we introduce 2 U(l) gauge fields, 
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and 5 M , and will consider effective Lagrangians with the following schematic structure in 
2+1 spacetime dimensions 



2 N v ., 

£ = { I C^A* — -^ci^)^^ 1 2 -|-s^ I | 2 J- """^Z^ { I C^a* — ^fe'At) ^« 1 2 -t-s^ | | 2 -I- ^^e^^^Ci^^^fo^-l- - ■ ■ , (1.1) 

a=l a=l ~ 

where //, u, A = x, y, r are spacetime indices, and s z and s v are the primary couplings we will 
tune to obtain our global phase diagrams. The integer k = 2 implements the needed semionic 
statistics. The ellipses represent additional terms in the effective potential for the z a and 
v a which are constrained by the projective symmetry group (PSG) i.e. the transformations 
of the spinons and visons under the symmetries of the lattice spin Hamiltonian. We will 
discuss these terms more carefully when we describe the different PSGs in the body of the 
paper. 

In passing, we note that supersymmetric versions of the doubled Chern-Simons the- 
ory in Eq. (11.11) have recently been the focus of intense interest in the string theory 
literature^^ 4 ^ 8 - 1 ^ 9 - 1 ^. Their model of interest is 4 ^ a Chern-Simons theory with a U(iV) 
x U(iV) gauge group, with opposite signs for the Chern-Simons term for the two U(iV)'s, 
and with matter fields which are bifundamentals in the U(iV)s. Eq. (11.11) is also precisely of 
this form with N — 1: we can define c M = a M + and = a M — and then the c M and 
dfj, fields have diagonal Chern-Simons terms with opposite signs, and the z a and v a carry 
bifundamental charges. The U(iV) x U(iV) theories have been argued 48,49 to be dual to M 
theory on AdS± x S 7 /Z) c , which is reason for the interest. The N = 1 case with M = 4 
supersymmetry has been argued^ to be exactly dual to a U(l) gauge theory without a 



Chern-Simons terms (the latter theory was reviewed in Ref. l5ll). The analogs of such N = 1 



dualities for non-supersymmetric theories are well-known in the condensed matter literature, 



and we will discuss examples in the present paper (see also Ref. l36l ). 

A natural question arises at this point: what are the conditions under which it is permis- 
sible to implement a U(l) CS theory realization of the Z2 spin liquid, rather than directly 
in terms of a Z 2 gauge theory? When we are discussing the topological properties of the 
ground state, or about single quasiparticle excitations, there does not appear to be any 
obstacle to using a U(l) theory 4 ^. However, the issue becomes more delicate when the 
excitations proliferate, and we are considering quantum phase transitions out of the spin 
liquid state. This question is discussed further in Section |Vll where we will find examples 
of transitions at which our U(l) CS description fails. However, we also find cases where it 
does succeed, and these are the main focus of this paper. As we will see in Section [TVt for 
these successful cases, because of the constraints of the lattice PSG, the lowest order terms 
which break either of the U(l) gauge invariances of Eq. (II. ip are of eighth order, ~ v\ as in 
Eq. (14.21) ; their effects are easily incorporated into our analysis as a soft symmetry breaking. 
We mention that the connection between doubled Z2 and U(l) CS theories has also been 
discussed by Balents and Fisher in a different context 52 . 
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Crucial to our analysis will be exact results on the low energy spectrum of C on a L x L 
torus as a function of s z and s v . For the case where both s v and s z are large and positive, 
both the spinons and visons are gapped, and we realize a Z 2 spin liquid. Here we can 
integrate out the z a and v a , and are left with a pure doubled Chern-Simons (CS) gauge 
theory. This theory was quantized exactly on a torus in Refs. I39IJ43I . The key variables in 



this quantization were the fluxes piercing the two cycles, C xy of the torus 



Ai= <p a^dx^ , B, t = <p b^dx^. (1.2) 
Jd JCi 

Given that all the matter fields carry unit or charges, the Ai and Bi should be regarded 
as periodic variables taking values on a circle of circumference 2n. After accounting for 
this periodicity, The solution of the ground state of the CS theory was found to be 4-fold 
degenerate; the degeneracy appears exponentially fast as L — > oo, provided the vison and 
spinon gaps remain finite. This 4-fold degeneracy is viewed as an essential characterization 
of the Z 2 spin liquicr^. 

The other phases in our phase diagrams appear when we allow one or both of s z and s v to 
vary to negative values. Then we can have phase transitions to new phases in which one or 
both of the z a and v a are "condensed" . However, the precise nature of the broken symmetry, 
if any, is not immediately obvious in such phases, given the presence of the 2 gauge fields and 
their CS term. The purpose of this paper is to describe these new phases and the associated 
quantum critical points. Here we note that the order parameter characterizing these states 
can be gleaned by carefully examining the low energy states of £ on a L x L torus. As an 
example, consider the state where s z is large and negative, and so a saddle point with z a ^ 
is favored. By global SU(2) spin symmetry, there are actually an infinite number of such 
saddle points along the manifold |^| 2 + 1^| 2 = constant i.e. along S 3 , the surface of a sphere 
in four dimensions. The low energy theory on a L x L torus can be expressed as a functional 
integral over S 3 , along with an integral over the gauge fields. We will solve this quantum 
theory exactly, and find an 'Anderson tower of states"—, with a non-degenerate ground 
state, and an infinite sequence of excited levels with energies ~ 1/L 2 above the ground state 
energy. The sequence of excited levels, and their degeneracies, can be uniquely identified 
with the quantum mechanics of a particle moving on S 3 /Z 2 , with the co-ordinates of the 
particle representing the average orientation of the order parameter across the entire torus. 
In other words, the primary effect of the gauge fluctuations is to reduce the order parameter 
characterizing the broken symmetry from S 3 to S 3 /Z 2 . It was these same gauge fluctuations 
which were responsible for the 4- fold degeneracy in the Z 2 spin liquid. The S 3 /Z 2 = 5*0(3) 
order parameter allows an immediate identification of the s z large and negative state: this 
is the spiral antiferromagnet, as found in CS2CUCI4. 

A related analysis for the other phases will be found in the body of the paper, allowing us 
to construct our phase diagrams. Here we show in Fig. [2] the phase diagram found for a case 
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FIG. 2: Global phase diagram for the case with N v = 1 and a particular model of the spinons 
(model Bill); a similar phase diagram appeared in Ref. @. The Neel antiferromagnet is found in 
a many materials, and has ordering wavevector Q = 27r(l,0). The geometry of the VBS state 
coincides with that foun d 20 ' 21 in EtMe3P[Pd(dmit)2]2- The spiral antiferromagnet shown in the 
figure has an ordering wavevector Q = 2tt(1 — e,0) with e = 1/6, small, as expected for J' < J. 
The experimentally realized spiral state in CS2CUCI4 has J' > J, and consequently a l arg er value of 
e « 1/2. The Z2 spin liquid in this phase diagram is similar to that proposed in Ref. |32j to explain 
observation o 26 ' 27 ' 28 ' 29 in k-(ET)2Cu2(CN)3. The transitions have been discussed previously: (2) 
the CP 1 field theory for the Neel- VBS transition^, (ii) the 0(4) field theory for the spiral- Z2 spin 
liquid transitio n 53 ' 54 ' 55 , (Hi) the mean field theory for the spiral-Neel transition^, and (iv) the 0(2) 
field theory for the VBS-Z2 spin liquid transition 7, ie . All these field theories are contained in our 
theory in Eq. (jl.ip . which also describes the multicritical point M. 



in which there is only one vison flavor, N v = 1 and for a particular model of the spinons — we 
label this theory model Bill. A mean-field phase diagram with the same phases appeared 
already in Ref. |8|, and here we shall show that these phases follow from the very general 
considerations outlined above, and also provide field theories for all the transitions and the 
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multicritical point M. It is encouraging that all the phases with broken symmetry, which 
descend from the Z 2 spin liquid with N v — 1, correspond precisely to those which have been 
experimentally observed so far. 

In Section [III we will describe the phase diagram of the theory in Eq. (11.11) as an abstract 
field theory, without reference to any underlying antiferromagnet. The specific spinon and 
vison degrees of freedom of the lattice antiferromagnet, and their possible PSGs, will be 
identified in Section IHII The combination of the results of Sections [II] and II I II lead to a 
variety of possible phase diagrams. These are described in Section HV] and the quantum 
phase transitions in Section |VJ Section IVII will give another semiclassical perspective on 
our results, which also identifies the limitations of the present U(l) CS approach. The 
concluding Section IVHI will make some further remarks on recent experiments. 

II. PHASES OF THE DOUBLED CHERN-SIMONS THEORY 

This section will discuss the phase diagram of the doubled CS theory, considered here 
as an abstract field theory. The interpretation of the phases in terms of the underlying 
antiferromagnet requires a more specific knowledge of the PSGs of the spinons and visons, 
and these will be considered in subsequent sections. 

For the case where the spinons and visons are gapped, as we have already noted, we 
obtain the Z 2 spin liquid. The fundamental property of this theory is the 4-fold degeneracy 
on a L x L torus, and this appears as a property of the pure doubled CS theory^ 1 ^. We 
are now interested in moving into one of the phases where one or both of the spinons and 
visons are "condensed" and understanding the nature of the broken symmetry. As in the 
Z 2 spin liquid, we will do this here by examining the low energy states of the theory on 
the L x L torus. We will compute the spectrum of Anderson "tower of states" with 
excitation energies ~ 1/L 2 : the spectrum of states will allow a unique identification of the 
ground state manifold (GSM) associated with the broken global symmetry. 

We will only consider here the case where there is a single spinon species, z, and a single 
vison species v; the generalization to the multiple species case is straightforward. For the 
broken symmetry phases, we need only consider the phases of these complex fields, and so 
we will write z ~ e %dz and v ~ e t6v where the corresponding symmetries are broken. 

As a warmup, consider first the case with a single broken U(l) symmetry, characterized 
by the U(l) order parameter e te , and no gauge fields. The low energy theory of 9 fluctuations 
is given by the action 



on a I x I torus, with 9 and 9 + 2tt identified. The couplings K 12 are two stiffnesses 
characterizing the broken symmetry. Because this is a Gaussian theory, we can make the 




(2.1) 
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mode expansion 



e(x,y,r) = e (r) 



2nmx 2-nny 



a k [r e 



ik-r 



(2.2) 



k^O 



where n and m are fixed integers (the winding numbers), 9$ represents the uniform fluctua- 
tion of the order parameter, and the a>k are the 'spin-wave' normal modes. Inserting this in 
the action we obtain 



S e = 47r 2 (m 2 + n 2 ) + / dr 



^(d T a k )* + fk 2 al 



(2.3) 



From this it is clear that the low-lying states have m = n = 0. The harmonic oscillators 
have energy ~ k ~ 1/L, while the 9q mode has energy ~ 1/L 2 . So for the lowest states, we 
put all the dfc oscillators in the ground state, and we obtain a tower of states with energy 



E„ 



E + 



2K X L 2 



(2.4) 



where p is an integer, measuring the angular momentum of the #o mode. 

It is now useful to note that the theory S$ in Eq. (12.1 j) is exactly dual to U(l) gauge theory 
with a Maxwell term, and so the latter should have the same tower of low energy states. 
Let us demonstrate this explicitly. First, we decouple the quadratic terms in Eq. (12.11) by 
an auxiliary current 



So 



d 2 rdr 



J 2 J 2 



2K X 2K 2 



(2.5) 



Integrating over 9, we obtain the constraint d^J^ = 0, which we solve by expressing J M in 
terms of a 'dual' gauge field a M : 

J„ = —t^xdydx- (2.6) 

Z7T 



The normalization of 1/ (2ir) is chosen so that periodicity of the flux variables A{ in Eq. (11. 2p 
with period 2tt is equivalent to the periodicity in the angular variable 9 — > 9 + 2ir. Now 
inserting Eq. (12.61) into (12.51) we obtain the U(l) gauge theory dual to Sg: 



S a 



d rdr 



8n 2 K 2 



%<K 2 K X 



(d x a y - d y a x f 



(2.7) 



where we have chosen the temporal gauge with a T = 0. An important consequence of the 
periodicity in variables is that the flux piercing the torus, J d 2 r(d. 



d y a x ) must be 



an integer multiple of 2tt. We can see this by moving the contour Cj in Eq. (II. 2p across 
the entire length of the torus: the change in the line integral upon returning to the initial 
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position must be an integer multiple of 2ir, and this change is equal by Stokes theorem to 
the flux piercing the torus. Thus the Hilbert space of S a breaks apart into distinct sectors 
with total flux 2rcp, where p is an integer. Within each sector, the ground state has zero 
photons (which are the dual of the spin waves), and has (d x a y — d y a x ) = (27ip)/L 2 . So the 
lowest energy state in each sector is 



Er, = E + 



2-rrp 
V ~T? 



E + 



2K l L 2 



(2.8) 



which is the same as the spectrum of Sg inEq. (12.41) . This verifies the equivalence of Eq. (12. lj) 
and (12771) . 

With these preliminaries out of the way, let us return to our doubled CS theory with one 
spinon and one vison. Consider the phase where the spinon is condensed and the vison is 
gapped, so s v ^> and s 2 < 0. Here we can simply integrate out the vison, and are left 
with the following low energy theory for 8 Z and the U(l) gauge fields 



S 7 



d rdr 



(d T e z 



K 



ik 



2tt 



(2.9) 



We can always choose the gauge 9 Z = (and A T = 0). In this gauge, the integral over a M is 
an ordinary Gaussian. Performing this integral, we obtain the action 



<S 2 



/ 



d rdr 



k 2 



8ir 2 K 2 



k 2 



8n 2 K- 



■(d x b y - d y b x f 



(2.10) 



Comparing this with the spectrum of S a in Eq. (12.81) . we obtain the low- lying states 



k 2 p 2 



E n = En 



2K\L 2 



(2.11) 



This shows that the theory S z in Eq. (12. 9p is equivalent to the U(l) scalar theory in Eq. (12. ip 
but with the periodicity 8 = 6 + 2ir/k. In other words, the GSM of this phase has been 
modified by the gauge fluctuations from S 1 to S 1 /Z k . Alternatively stated, the broken 
symmetry of the ground state is associated with distinct values of the composite field z k . 

It is useful to have another perspective on the above result by an alternative analysis of 
the theory S z in Eq. (12. 9p . This analysis begins by 'undualizing' the gauge field into a 
dual scalar 6b- For this we introduce, as in Eq. (12.61) . the current = e^\d u b\/ \2n) and 
impose the constraint <9 M = by a Lagrange multiplier 6^, this modifies S z in Eq. (12.91) to 



d 2 rdr 



K K J b2 

-^(d T 6 z - a T ) 2 + -£(dA + a t ) 2 + U^ka, - d fl 6 b ) + 
* ^ 2K 



[212) 



The last term is a useful regularization, and the original theory in Eq. (I2.9p is obtained in 
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FIG. 3: Schematic phase diagram of the doubled CS theory in Eq. (jl.ip for k = 2. All non-zero 
order parameters which characterize the broken symmetry in each phase are shown. 



the limit K — > oo. Now we perform the integral over J b and obtain the theory 



S 7 = / d rdr 



^-(d T 9 z - a T f + ^{dA + a t ) 2 + ^{ka, - d,9 b f 



(2.13) 



In the limit K — ► oo, we see that we must have a M = (l/k)dp9i,. However, 9b is a variable 
periodic under 9b ^ 9 b + 2n, and hence the periodic flux variables Aj in Eq. (jl.2p can only 
take the values 

A = 2 -f (2.14) 

where the pj are integers. In other words, U(l) gauge field a M has been reduced to a Z k 
gauge field. Thus e lke ~ z fc is gauge invariant, and this explains our results above on the 
distinct values of z k identifying distinct ground states. 

We have now completed our discussion of the state where the spinon is condensed and 
the vison is gapped (s v ^> and s z C 0). Clearly, the complementary phase where the 
vison is condensed and the spinon is gapped (s v <C and s z 3> 0) is amenable to a parallel 
treatment, with complementary results and a v k order parameter. 

Finally, let us consider the case where both the visons and spinons are condensed, s v <C 
and s z <C 0. In this case, we can see in the gauge 9 Z = and 9 V = that both fields a M and 
6 M are fully gapped. So there is a unique ground state, and no other excited states whose 
energy vanishes as L —>■ oo. 

We can now generalize these results to the cases with multiple flavors of visons and 
spinons, and the results are summarized in Fig. [3] for k = 2. The flavor indices simply tag 
along for the order parameters involving a k-fo\d composites of spinons or visons, which 
are invariant under the gauge transformations. However, there are now also additional 
gauge-neutral order parameters possible, like z^zp, which were absent for the single flavor 
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case, because they were not associated with any broken symmetry. 

Armed with the results in Fig. [31 and with a knowledge of the microscopic PSGs of 
the spinons and visons (which are described next in Section UTTj) . we can easily deduce the 
physical characteristics of the phases of a variety of antiferromagnets. 

III. SPINONS AND VISONS 

A rich variety of spinon and vison operators can be defined for S = 1/2 antiferromagnets 
on the lattice in Fig. HJ and we shall not attempt any complete classification. Clearly, 
the choice depends sensitively on the details of the microscopic Hamiltonian. However, in 
previous semiclassical analyses, a few natural choices have emerged for different limiting 
values of J' / J. We will describe these below, and show how they can fit together in a 
doubled CS theory like Eq. fTO]) . 

The essential characteristics of the spinons and visons will be their transformations under 
the symmetry operations of the underlying spin model. These symmetries are the lattice 
translations T\ and T2, the lattice reflections P x and P y , and time reversal, T: 



Ti 


(x 


V)- 


(x+ l,y) 


T 2 


(x 


V)- 


(x+l/2,y + VS/2) 


P x 


(x 


I/)- 




Py 


(x 


V)- 


(x, -y) 


T 


t - 


-> -t 





Spin rotation is also a symmetry, and is easily implemented by contracting the spinor indices. 
The following subsections will consider the PSGs of spinons and visons in turn. 

A. Spinons 

One natural model of spinons appears upon describing the spiral ground state of the tri- 
angular lattice antiferromagnet in terms of Schwinger bosons. So we write the spin operator 
on the lattice sites as S = b'ab/2, where a are the Pauli matrices. For the Schwinger bosons 
we make the following low energy expansion 57 in terms of the spinon fields, z a : 

b a ~ z a exp(iQ ■ r/2) + ie a/9 ^exp(-zQ • r/2), (3.2) 

From this parameterization we can then deduce the following expression for the spin oper- 
ators 

S — Hi cos(Q • r) + n 2 sin(Q ■ r), 
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z — 


> e lQ ^ 2 z 




z — 






z a 


* 




z — 




T 


z — 


> zz*. a 



ni = Re[z a v az], ti 2 = lm[z a v az\, 

n 3 = Hi x n 2 = z^az. (3.3) 

We observe that Q is the ordering wavevector of the spiral. So for J' « J, we expect 
Q « (27r/3, 0). For CS2CUCI4, which has J' ~ 3 J, we have Q ~ (it, 0). Finally, in the square 
lattice limit, J' <C J, we have Q ~ (2tt, 0). Thus we expect Q to increase monotonically from 
(it, 0) to (2tt, 0) with decreasing J'/ J. Also note that 711,2,3 are three mutually orthogonal 
vectors. 

The parameterization in Eq. (13.31) allows us to deduce the PSG of the z a . These z a 
spinons will couple minimally to the a M gauge field, and so there is a natural implied PSG 
for the a M . We call the resulting PSG of spinons as model A: 

Spinons, Mo del A (3.4) 

<V -> <V 

tin 

a x — > a x , a y — > —a y , at — ► ~ a t 
> a x , a y — > — a y , — ► 

-> a M . (3.5) 

We note that under this model A PSG, ti\^ are odd under time reversal, while ti% is even. 
From this, and the representation in Eq. ( 13. 31) . we deduce that n3 is a spin nematic order 
parameter. 

The Model A spinons are natural for J' > J, where the spiral order is likely to be 
present. However as we approach the square lattice limit with J' — > 0, there is a finite 
range of small J' < J over which we expect that Q is pinned exactly at (2tt, 0), and we 
have the conventional 2 sublattice Neel order appropriate for the square lattice. In this 
case sin(Q • r) = identically, and cos(Q • r) = cos(27nc) = (— l) 2x . In this limit, we 
can define another model of spinons which appeared in previous theories of square lattice 
antiferromagnetsi^. We map z a — ► (z a + ie a pz%) /a/2 and then find that Eq. (13.31) is replaced 
by 

S ? = m 1 (-l) 2x , mi = z^az, (3.6) 

A significant difference from Eq. (13.31) is that now the U(l) gauge invariance associated with 
a M is explicit, because the representation (13.61) is invariant under the gauge transformation 
z a — > z a e ld . We label these spinons model B, and they also have mappings under the square 
lattice space group, which we can deduce from Eq. (13. 6p to be: 

Spinons, Model B (3.7) 
Ti : z -> -z, a M -> 
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T 2 

p x 

Py 
T 



z — > iz, a x — 

Z > Z, d x * 



x > Q"y 



a 



► -a. 



yi ®t 

yi 



a t 



(3.8) 



It is now clear that under the model B PSG, rh\ is odd under time-reversal. In the state 
where the spinons are condensed and the visons are gapped, we see from Fig. [3] that we also 
have two additional vectors in spin space which characterize the broken symmetry in the 
ground state (analogous to the 3 vectors found in model A): 



rh 2 + irriz 



z t a v az 



(3.9) 



These vectors do not appear in the present expression for the spin operator in Eq. (13.61) . 
and so their physical interpretation is not yet clear. Let us, therefore, compute the PSG of 
these vectors: 



Model B 



(3.10) 



Tl 

p x 

T 



m 2 ,3 
m 2 - 

m 2 - 



"> "^2,3 
m 2 77T.3 

-> -m 2 ,3 

-* ^2,3 

fh 2 m,3 



-m 3 



-m 3 



(3.11) 



From these properties it is evident that we can identify m 2 as the central axis (in spin 
space) about which the spins are precessing in the spiral antiferromagnet shown in Fig. [21 
We will present a more detailed analysis in Section IIVDI which shows how the spiral state 
emerges for model B spinons. This identification is also consistent with the analysis of this 
phase in Ref. @, where the spiral phase was induced by the condensation of a charge 2 Higgs 
scalar — the vison is dual to this scalar, and the gapping of the vison is equivalent to the 
condensation of the Higgs scalar. 



B. Visons 

In the simplest models of visons^k^, we take real particles hopping on the sites of the 
dual lattice, subject to a flux of n around every site of the direct lattice. In other words, the 
vison is the Ising field of an Ising model on the dual lattice, with exchange couplings chosen 
so that every plaquette surrounding a direct lattice site is frustrated. For the antiferromagnet 
on the lattice in Fig. [TJ the Ising model resides on the dual lattice shown in Fig. HI Because 
of the dual relation between the couplings of the antiferromagnet, and the model of the 
visons, we expect that w'/ w decreases as J'/ J increases. 
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FIG. 4: The dual Ising model describing the vison dynamics in the dual honeycomb lattice. All the 
vertical bonds have hopping amplitude to', all the other bonds have hopping amplitude w. Note 
that a small J'/ J implies a large w'/w, and vice versa. Besides the lattice anisotropy, there is 
a 7T— flux through every hexagon, which frustrates the vison kinetics; this flux is implemented by 
changing the sign of the dotted w bonds. 



It is a relatively straightforward matter to obtain the spectrum of such a particle moving 
on the lattice in Fig. HI Because the kinetics of the vison is frustrated by the background 
spinon charge on every site, the product of the hopping amplitudes on the six links around 
each hexagon is -1. We are free to choose one of the links to be negative, and our convention is 
shown in Fig.Hl The spectrum has multiple minima in the Brillouin zone, and we introduce a 
real vison field for each such minimum in the spectrum. This procedure parallels that carried 



out in obtaining the multiple vortex flavors in Refs. |58| and |59|, but with modification that 
we are considering real vison fields and not complex vortex fields. 

For the general set of parameters for the lattice in Fig. HI we find that there are either 
4 or 2 minima of the vison dispersion lattice in the Brillouin zone. The 4 minima occur 
near J' « J, and so are appropriate for the triangular lattice limit; indeed for J' = J these 
minima co-incide with those found by Moessner and Sondhi^. For J'/ J small (w'/w large), 
near the square lattice limit, we find only 2 minima. 

Let us begin by considering the 4 minima case. These are at momenta of the form 

qi = (gix,vr/2), 

q 2 = (tt - q lx ,n/2), 

q3 = (-^,-71-/2), 

q 4 = (q lx -7r,-n/2). (3.12) 
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There are two choices to combine these four real minima to two complex minima, which will 
then correspond to the complex vison fields v a , with a = 1,2. The first choice, which we call 
model I, is: 



vi ■ qi = (gix,7r/2), 

v 2 : q 2 = (vr - q lx ,n/2), 

K- qs = (-qix,-^/2), 

v* 2 ■ q4 = (gi x -7T,-7r/2), (3.13) 



while model II is: 



«i : qi = (Qix,n/2), 

V*2 ■ q2 = (tt - gia.,7r/2), 

: q3 = (-Six, -t/2), 

^2 : q4 = (qix - 7T, -7r/2). (3.14) 



These two choices lead to 2 models for the vison PSGs: 



Visons, Model I 





Ul - 




T 2 


Vi - 




P. 


vi - 


-> 




vi - 




T 


V a - 





v 2 -> -e l7 t>2, frr -> 6 X , ^ -> -63,, 6* -> -6 t 
-> uj, b x -> -6 X , &„ -> &„, 6t -> -6t 

(3.15) 



where 7 and are two incommensurate angles. 
The second model for the vison PSG is 



by, b t -> -b t 
b t 

(3.16) 

where again 7 and are two incommensurate angles when the lattice is distorted. For an 
isotropic triangular lattice, 6 = —tt/12, 7 = ir/6, and they will continuously evolve to vison 
model III by tuning the distortion of the lattice. 

Finally, let us move to the case where the four minima in the vison band merge to two. 
In model I, v\ and v 2 merge together, while in model II v\ and v\ merge together. But v\ 
and v 2 carry opposite gauge charges in the CS theory, and so this merger violates the gauge 



Visons, Model II 
Tl . ~ ~ e m* Vu V2 _ e -™+*n* V2i b^b^ 

T 2 : vi -> e ie v 2 , v 2 -> ie i6 v u & M -> 6 M , 

P x vi —> e-^vl v 2 - -e'^v*, b x - 6,, &„ - - 

P y : fi -> v 2 , v 2 —>vi, b x ->b x , b y -b y , b t -> 

T : v a -> v* a , &„ -> &„. 
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invariance. So if we want to evolve smoothly from four minima to two minima, we have to 
take model I. 

After the merger, the two minima are located at (7r/2,7r/2) and (— n/2, — 7r/2). So we 
have just to use one component complex vison v, and its PSG leads to 

Visons, Model III 
7i : v -> iv, 6 M -> 
T 2 : t; - -e^/S*, &„ - -&„ 
P x : f -> -zV, -> 6 y -> -by, b t -> -6 t 
Py : u -> fer -> -&x, ^ — > fey, b t -> -fet 

T : w ^ w*, (3.17) 

Finally, we consider the nature of the vison order parameters v a Vb and v*Vb which appear 
in the phases in Fig. [3J 

The simplest case is model III, with only one complex vison v, in which case the only 
non-trivial order parameter is v 2 . From the PSG, we see that v 2 is the square lattice VBS 
order parameter, associated with the VBS state shown in Fig. [2j Using the definitions Vx,y 
for this order parameter in Ref. 0, we have Vx ~ cos(0 + vr/4), Vy ~ cos(0 — tt/4) where 
v 2 = exp(z0). Here (and henceforth), the axes x, y refer to the principle axes of the "square" 
lattice formed by the J bonds in Fig. [lj 

The vison order parameters for the other models also describe VBS orders but of a differ- 
ent nature. The vison operator v a is subject to a Z 2 gauge invariance, therefore the physical 
VBS order parameter should always be bilinear of v a . There are in total 15 independent 
bilinear of v a , and the detailed VBS pattern drive by vison proliferation depends on the 
Hamiltonian of visons, which will be discussed in the next section. 



IV. PHASE DIAGRAMS 

Now we turn to the crucial question of combining the spinon and vison PSGs in Section [1111 
into consistent theories of the form in Eq. (11.11) . We will denote the resulting theories by an 
obvious notation i.e. the theory Bill has spinons under model B and visons under model 
III. 

In principle, there are now 6 possible theories, AI, All, AIII, BI, BII, and Bill, and 
associated phase diagrams. To establish the consistency of these theories, we have to examine 
the transformation of the CS term under the respective spinon and vison PSGs. The results 
of such an analysis are summarized in Table [I] for all theories. We see that under theories 
BI and Bill the Chern Simons term is strictly invariant, and so these theories are clearly 
consistent. For the remaining theories, the overall form of the CS term remains invariant, 
but some of the transformations do lead to a change in sign of the CS term. However, the 
role of the CS term here for k = 2 is only to implement a mutual semionic phase of tt, 
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BI and Bill 


AI and AIII 


All 


BII 


J 1 


+ 


+ 


+ 


+ 


T 2 


+ 




+ 




P x 


+ 






+ 


Py 


+ 


+ 






T 


+ 






+ 



TABLE I: Transformation of the mutual Chern Simons term under the space group operations for 
the various theories. The CS term changes its overall sign, as indicated. 

and this is invariant under the sign change. Equivalently, we are free to define the vison 
at momentum Q to be either v or v*, which means that in the system there should be 
particle-hole symmetry of vison i.e. on average the spinons see zero flux. This particle-hole 
symmetry corresponds to the free choice of the sign of gauge charge of vison, and leads to 
the freedom of the sign of the mutual CS term. 



A. Model AI 

The Lagrangian should be invariant under all the symmetry and PSG transformations, 
which in general takes the form: 

2 2 ., 



a=l a=l 

2 



+ U z \z a \\ + Uv \Vg\\ 



+ g\vi\ 2 \v 2 \ 2 + ••• (4.1) 



Let us first identify all the symmetries of this Lagrangian. The U(l) gauge symmetries 
associated with gauge field a M and b^, correspond to two global U(l) symmetries U(l) a x 
U(l)b in the dual picture, which lead to the conservation of gauge fluxes. Through the 
mutual CS term, the gauge flux of a M is attached with the vison number, and the gauge 
flux of b^ is attached with the spinon number. On top of the U(l) gauge symmetries, the 
global symmetry of this Lagrangian to the fourth order of z a and v a is SU(2) spin x U(l)xZ 2 . 
The U(l) symmetry corresponds to the U(l) transformation on vison bilinear v\v2\ the Z 2 
symmetry corresponds to interchanging V\ and v 2 , which physically can be understood as 
the reflection symmetry P y . If the lattice is an undistorted triangular lattice, g = 0, and 
the vison doublet enjoys an enlarged SU(2) flavor symmetry in this mutual CS theory (or 
an 0(4) symmetry in a theory with only vison). The ellipses in Eq. (14.1 1) include terms no 
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FIG. 5: Phase diagram of model AI with g > 0. The stripe order is illustrated in Fig. [7J while the 
spiral antiferromagnet is as in Fig. [2j 
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FIG. 6: Phase diagram of model AI with g < 0. The square lattice VBS state is as illustrated in 
Fig. El 



less than sixth order of v a , which may introduce higher order anisotropy. In the distorted 
triangular lattice, the lowest order term which breaks the symmetries in Eq. (14.11) is at the 
eighth order: 

£ 8 = <78(^2) 4 + H.c. (4.2) 

In the undistorted lattice, the lowest order symmetry breaking term is at the sixth order. 
If only the terms below fourth order are considered, we can minimize the Lagrangian in 
Eq. (14. ip with tuning parameter s z and s v , and obtain the phase diagrams in Fig. O and El 
with the phases described in the following subsections. 



1. The Z2 spin liquid and spiral phases 

The phase with both the spinons and visons gapped out (s z > 0, s v > 0) is the Z 2 
spin liquid, as was discussed in Sections fl] and [TTJ has four-fold topological degeneracy on 
a compact torus. The phase with visons gapped and spinons condensed (s v > 0, s z < 0), 
is the incommensurate spin spiral state, with wave vector Q. With the vison gapped, the 
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FIG. 7: Illustration of the Z 2 stripe order. This can be viewed as a bond density wave, with 
alternating signs on successive rows of bonds. 



gauge field 6 M is in the photon phase, and so the CS term "Higgses" out the gauge field a M . 
Or more precisely, for k = 2, one gauge flux of carries two gauge charges of a M ; therefore 
the photon phase of which is the superffuid phase of gauge flux, breaks the U(l) gauge 
invariance of a M to Z 2 , as was discussed in Section [Til As also noted in Section [TTl this implies 
that the GSM of the spinon condensate is S 3 /Z 2 , and this is the GSM of the spiral spin 
state. Another way to understand this phase is that, because the vison number is attached 
with the flux number of a M through the mutual CS term, in the Mott insulator phase of 
vison the photon of a M is gapped out, while is in the photon phase. The global symmetry 
U(l)b becomes the global U(l) symmetry of the spinon z a , which according to the PSG in 
Eq. (13.51) corresponds to the physical translation transformation. 

2. The vison condensate with s v < and s z > 

The nature of the phase with spinon gapped and vison condensed depends on the sign 
of g in Eq. ( 14.11) . With the spinon gapped, the U(l) gauge field b^ is broken down to Z 2 
gauge field. Integrating out the remnant Z 2 gauge field, the vison v a enjoys a U(l) x U(l) x 
Z 2 symmetry. The two U(l)s correspond to the global symmetry of two flavors of visons 
respectively, and the Z 2 symmetry corresponds to the interchange symmetry between V\ 
and v 2 . With g > 0, the vison condensate breaks the U(l) x U(l) x Z 2 symmetry to U(l) 
symmetry i.e. only one flavor of v a condenses. Let us assume z\ condenses, and z 2 remains 
gapped. The GSM of the vison condensate is S 1 x Z 2 . The Z 2 degeneracy is described by the 
Ising order parameter v ^a z v , which corresponds to the stripe order depicted in Fig. [71 The S 1 
in the GSM is described by the order parameter vf, which corresponds to an incommensurate 
valence bond density wave along the x direction with wave vector (2g la .,7r), meanwhile 
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breaks the Ising symmetry of interchanging V\ and V2, or the reflection symmetry P y . The 
continuous symmetry U(l) transformation of V\ which has been broken by this valence bond 
density wave is the translation along x. Because V\ carries an incommensurate momentum, 
the full Lagrangian in (14. ip with all the higher order perturbation should preserve this global 
U(l) symmetry, and the GSM S 1 x Z2 is not broken down to smaller manifolds. For instance, 
the eighth order term £ 8 violates the vison numbers of both v% and t>2, in the phase with 
only v\ condensed, £ 8 is suppressed. 

If g < 0, the vison condensate breaks its global symmetry to Z 2 i.e. both v\ and t>2 
condense with equal stiffness, the GSM is S 1 x S 1 if there are no more symmetry breaking 
terms. The two S 1 corresponds to two U(l) order parameters V1V2 and v\v2 respectively. 
The U(l) symmetry of v\v 2 is preserved by the full Lagrangian, while the U(l) symmetry 
of V\v 2 is broken by the eighth order term £ 8 i n Eq. (j4.2|) . This term breaks the U(l) 
symmetry of v \v 2 to Z4 symmetry. Actually, the order parameter V1V2 corresponds to the 
degeneracy of the four VBS states, which are smoothly connected to the four fold degenerate 
VBS states in the square lattice limit with J 1 ~ 0. Using the square lattice coordinates, the 
VBS order parameters are 

Vx ~ V\V 2 exp(i7r/4) + v*v% exp(— wr/4), 

Vy ~ v 1V2 exp(— in/A) + t>*t>2 exp(i7r/4). (4.3) 

The square lattice VBS order is selected when g 8 > 0, otherwise the four fold degenerate 
plaquette state is favored. On top of this commensurate VBS order, an incommensurate 
valence bond density wave corresponding to v\v 2 is also present, with wave vector (n — 
2q lx ,0). 

3. The phase with both spinous and visons condensed, s z < 0, s v < 

The most interesting phase is the phase with both spinons and visons condensed. In this 
phase, the physical order parameter should be the U(l) gauge invariant bilinears of z a and 
v a , as discussed in Section [TXJ The spin order parameter is the nematic vector ~ z^a a z. 
The VBS pattern, depending on the sign of g, is either the Z2 order parameter v^o~ z v , or 
the incommensurate valence bond density wave v*V2- Note however, that the commensurate 
VBS order parameters Vx and Vy vanish, because they are not gauge invariant. Another 
way of understanding the vanishing of V^ and Vy is as follows: when the spinon z a is still 
condensed, the flux of a M is in the Mott insulator phase; because the flux number of a M is 
attached to the vison number through the mutual CS term, any order parameter violating 
the vison number conservation should not condense. 

The GSM of this phase is either S 2 spin x Z 2 (g > 0) or S 2 spm x S 1 (g < 0). And these GSM 
are not lifted by any higher order term in the Lagrangian Eq. 14.11 

Also note that £ 8 in Eq. (14.21) violates the gauge symmetry of b^. In the condensate 
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FIG. 8: Phase diagram of model All with g < 0. The phase diagram of model All with g > is 
the same as that for model AI with g > 0. 



of visons, Lg confines the fluxes of 6^, and hence the spinon excitation z a is also confined, 
which is consistent with the intuitive understanding of VBS states. 

B. Model All 

The phase diagram of this model is very similar to the previous subsection, the model 
AI. The only difference is that we now replace v 2 by v 2 . So the spiral spin density wave, the 
Z 2 spin liquid, and the VBS order are the same as in model AI. The commensurate VBS 
order parameter is now represented as 

Vx ~ V\V 2 exp(i7r/4) + v\v 2 exp(— wr/4), 

Vy ~ Viv'2 exp(— in/A) + v\v 2 exp(i7r/4). (4.4) 

And the incommensurate valence bond wave is represented by v \V 2 . 

In model All, the VBS pattern with both spinons and visons condensed is different from 
model AI. Because in this phase all the physical order parameters should be U(l) gauge 
invariant, the VBS order parameter is either the Z 2 symmetry breaking v'a z v, or the Z 4 
symmetry breaking V% and Vy depending on the sign of g. The GSM is S 2 pin x Z 2 (g > 0) 
or Sg pin x Z4 (g < 0). The phase diagram is shown in Fig. [HJ 

C. Model AIII 

In this model there is only one complex vison. The Z 2 spin liquid and the spiral spin 
state are the same as the two previous models. The vison condensate induces the four fold 
degenerate VBS order, with order parameter V^ ~ v 2 exp(i7r/4) + h.c, Vy ~ v 2 exp(— m/4) + 
h.c. However, one can no longer write down a U(l) gauge invariant order parameter in terms 
of v, therefore the phase with both spinons and visons condensed has only the nematic order 
r?3, and no other lattice symmetry breaking. The phase diagram is shown in Fig. [91 
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FIG. 9: Phase diagram of model AIII. 



If the parameter t'/t is tuned, the four vison minima will merge to two vison minima i.e. 
the two complex visons become one complex vison. Therefore by tuning t'/t, model AIII 
can be connected to model AI. 

D. Models BI, BII, Bill 

These models are similar to models AI, All, AIII. The main difference is that, in the 
phase with both spinons and visons condensed, the spin order is the collinear Neel order, in 
place of the spin nematic order parameter. The phase diagram for model Bill was shown in 



We now discuss in some detail how the spiral order emerges in model B, as this is not 
evident from the underlying spin representation in Eq. (13. 6p . If we use only the constraints 
imposed by the model B PSG in Eq. (13. 8p . then the Lagrangian of the spinons allows an 
additional inear spatial derivative term: In the Lagrangian of spinon of 



The term C x violates the enlarged U(l) gauge invariance of the mutual CS Lagrangian 
discussed below Eq. (13. 6j) . The mutual CS term will bind this term with the monopole 
operator of 6 M , which creates 2n gauge flux. We denote this monopole operator by A4 b , 
then in the U(l) gauge invariant formalism C x reads 



In the phase with both spinons and visons condensed, the term C x is suppressed, because 
of the conservation of the flux of which is attached to the spinon number of z a . However, 
once the visons are gapped, the monopoles M. b condense and C x becomes relevant. Due to 
its linear derivative of x, C x will drive the system into an incommensurate spiral state with 
wave vector along x axis, as has been described in Ref. @. The size of the incommensurate 



Fig. El 



X 



e a pz a V x zp + H.c. 



(4.5) 



C x ~ M b e a pz a V x zp + H.c. 



(4.6) 
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wave vector increases linearly with (Mb) ~ y/P^ wnere Pb is t ne stiffness of the flux 
condensate, which is proportional to the gap of vison. 



V. QUANTUM PHASE TRANSITIONS 

There are many phase transitions involved in the phase diagrams discussed in Section | 
We will study them in the same manner as the previous section. 

Before turning to the individual cases, it is useful to discuss an alternative form of the 
mutual CS theories, Eq. (11.11) . For many of the vison models, it is possible to^>^ "undualize" 
the vison degrees of freedom: this leads to an alternative formulation of the theory, now 
without a CS term. This new undualized form will be useful for many purposes. 

Let us first consider the simplest case of a single complex vison, as in model III. By the 
usual boson-vortex duality^, the dual of v is the monopole operator Mb introduced below 
Eq. (14.61) . This monopole operator carries charge k = 2 under a M , and consequently we can 
write the theory for the two-component spinor z a and the complex "Higgs" scalar Mb- Thus 
a theory equivalent to Eq. (11.11) for models AIII and Bill is 



z 

= ^2{\(d^-ia^)z a \ 2 + s z \z a \ 2 ^ + \(dp + 2ia^)M b \ 2 - s v \M 

a=l 

+ u z i^2\z a \\ +u M \M b \ 4 + v M \M b \ 2 



|2 



2 

|2 



,a=l 



+ \{M b e a pZ a d x zp + H.c.). (5.1) 

The A term descends from Eq. (14.61) . and is present only for model Bill. A closely related 
model, for a similar model, was obtained directly from the Schwinger boson formulation in 
Ref.fl Note that we have (schematically) changed the sign of the "mass" term for M b from 
that for the vison v. This reflects the dual relation between the fields, and the fact the v 
is condensed when is gapped, and vice versa. We note that the mapping between the 
CS theory in Eq. (II. ip and the non-CS theory in Eq. (15. ip is similar to that described for 
super symmetric gauge theories in Ref. 0. 

A similar (un) duality mapping can applied to visons in model I and II. This mapping 
only works for the g < ("easy-plane") case of the theory in Eq. (14.11) . In this case, the 
vison fields v y and the gauge field form an easy-plane CP 1 model, and so we can directly 



use the duality mappings of Ref. |6l|. The dual theory is yet another CP 1 model, with fields 
mi and irt2 and a gauge field c^. Here my are merons in the vison CP 1 model, and the 
monopole in the b^ field is^i Ai b ~ m\m2- Thus a new form of the theory (14.11) for models 
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AI, All, BI, BII with g < is 



z 

a=l 

+ + ia^ + ic^rail 2 + I (<9 M + ia^ - zc At )m 2 | 2 - ^ j|mi| 2 + |m 2 | 2 j 

+ u z (^y]\Za\^ +M m ||mi| 2 + |m 2 | 2 J +# m |mi| 2 |m 2 | 2 

+ w m ||mi| 2 + |m 2 | 2 | ^2 \ z a\ 2 ^j + X(m 1 m 2 e al3 z a d x z f3 + H.c). (5.2) 

Again the A term is present only for models BI and BII. Also the phase diagrams can be 
mapped by keeping in mind the dual relation between the visons t>i j2 and the merons m 12 : 
the visons are condensed when the merons are gapped, and vice versa. 

Now we will turn to a description of the transitions for the various models, using the 
theories in Eq. ([HI]), floTTD and (fOft . 

A. Phase transitions in model AI 

1. Transition between Z 2 spin liquid and spiral spin state 

This transition is known^^^ to be a 3D 0(4) transition, and the mutual CS theory does 
reproduce this 0(4) universality class: in the Z 2 spin liquid, the vison is gapped, therefore 
the gauge field is "Higgsed" by gauge field 6 M , and the U(l) gauge symmetry of a M is 
broken down to the Z 2 gauge symmetry. The critical point described by spinon z a enjoys 
an enlarged 0(4) symmetry. 



2. Transition between Z 2 spin liquid and the VBS state 



The nature of this transition depends on the sign of g. When g < 0, the Lagrangian 
describing this transition is 

2 

£ = ^ \ d H V »\* + r \ V a? + Uv\Va\ 4 + (2u v + </) \v 1 1 2 \v 2 1 2 + </ 8 «l«2- ( 5 ' 3 ) 
a=l 

This Lagrangian describes two coupled 3D XY transitions. The coupling g$ is clearly irrel- 
evant at this 3D XY transition. The scaling dimension of u\ 2 = 2u v + g is 2/u — D < 0, 
therefore is also irrelevant [y is the critical exponent defined as £ ~ r~ v at the 3D XY 
transition, which is greater than 2/3). So the transition between the Z 2 spin liquid and the 
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VBS order is two copies of 3D XY transitions when g < 0. 

When g > 0, the transition breaks the U(l) x Z 2 symmetry. There can be one single first 
order transition or two separate transitions, with 3D XY and 3D Ising universality class 
respectively. If the triangular lattice is undistorted, g = 0, there is one single transition 
between the Z 2 spin liquid and the VBS order, which belongs to 3D 0(4) universality class. 

3. Transition between spiral spin state and the nematic+VBS state 

If g < (g > 0), this transition is described by a CP 1 Lagrangian with easy plane (easy 
axis) limit: 



The eighth order anisotropy term £ 8 is suppressed at this transition, because the condensate 
of spinon is the Mott insulator phase of the flux of a M , which guarantees the conservation of 
total vison number. 

4- Transition between VBS order and nematic order 

This is a CP 1 transition described by spinon z a and U(l) gauge field a M . The eighth 
order anisotropy term of vison v\v\ + h.c. violates the conservation of the flux of a M i.e. it 
corresponds to the instantons in the 2+ld space time which creates/anihilates gauge fluxes. 
Because k = 2 in Eq. (14. ip . one flux of a M carries two visons, therefore £ 8 corresponds to 
a quadrupole process. If g > 0, only one component of V\ and t> 2 condenses, the quadrupole 
process which involves both v\ and v 2 is suppressed. However, if g < 0, the quadrupole pro- 
cess is present, but expected to be irrelevant at the CP 1 critical pointA We note here recent 
numerical studies of the CP 1 field theories, which include indications that this transition is 
weakly first ordepSLSlS&S&SML. 

One other issue to notice is that the spinon velocity and the vison velocity do not have 
to be equal. Therefore in the CP 1 models described above, the velocity of matter fields and 
the velocity of gauge fields are essentially different. In the large N limit, the U(l) gauge 
field has scaling dimension 1, and the one- loop self-energy of gauge field leads to the same 
velocity as the matter fields. The term with velocity anisotropy has scaling dimension 4, 
and hence is irrelevant for large enough N. 



2 




(5.4) 



a=l 
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B. Phase transitions in model All 

The phase transitions in model All are similar to model AI. The only difference is that 
the eighth order vison term C 8 = g^iyiV^Y + h.c. conserves the total flux number of a M , 
therefore there is no quadrupole process at the transition between the nematic order and 
the VBS order. 

C. Phase transitions in model AIII 

In model AIII, the transition between the spiral spin order and the Z2 spin liquid is 
still 0(4), while the transition between the Z 2 spin liquid and the VBS order is a 3D XY 
transition (with an irrelevant eighth order anisotropy), because there is only one flavor of 
vison. 

The transition between the spiral spin order and the nematic spin order is an inverted 
3D XY transition 60 , or a CP model with one component of complex boson v coupled with 
the U(l) gauge field b^. The transition between the nematic order and the VBS order is a 
CP 1 transition with irrelevant quadrupoles. 

D. Phase transitions in model BI, BII and Bill 

The transitions in the models with spinon B are similar to the models with spinon A. 

The major concern is the effect of C x in Eq. (14. 6p at the critical points. As discussed in 

Section IIVD[ this term is only effective when the vison is gapped or critical, for instance 

the transition between the Z2 spin liquid and the spiral antiferromagnet. The theory for 

this transition with spinon A was the 0(4) model. With the term C x present in model B, 

we can redefine the spinon field using a x-dependent 0(4) rotation to absorb the linear x 

derivative^. The transition is therefore seen to remain in the 0(4) class in model B. 

At the transition between the Neel order and spiral order, the field theory is given by a 
C p(iv-i) model with N flavors of 

visons (N = 1 or 2) and gauge field 6 M . C x violates the 
conservation of spinon, and hence corresponds to a monopole term A4b of b^. For simplicity, 
let us consider the model Bill with one vison component as an example. Here we can analyze 
the Neel-spiral transition from the theory in Eq. (15.11) by condensing the monopole operator 
A4b- We parametrize the spinon z as z = e ia {e 1 ^' 2 cos(0/2), e~ 1 ^' 2 sin(6 l /2))', where a is a 
gauge dependent phase angle coupled with the gauge field a M . Then the effective Lagrangian 
can be written as 

L = {d^f + (sinfl) 2 (,9 M 0) 2 + ~\(V x 9Re[M b ) + V x <f>emffhn[M b ]). (5.5) 

where M.^ is the gauge invariant monopole M.\> = M.b^ 2ia - Integrating out the gap less 
spin-waves and 9, a singular long range dipole interaction is generated for field M. b with 
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momentum dependence q 2 /{q 2 + uj 2 ), which will change the relative scaling dimension be- 
tween x and y, r. The effective theory for XY field ^ ~ Aib can be viewed as an effective 
z = 2 theory with scaling dimension A[q x ] = 2A[q y ] = 2 A [a;] = 2: 

u = -rt~ 2 m 2 + (q 2 y + u 2 )m 2 + 9 m 4 + ■ ■ • (5.6) 

qy + uj 2 y 

The upper critical dimension of this z = 2 field theory is d — 2, therefore this transition will 
be a mean field transition instead of a 3D XY transition. 

For N v = 2 theories in Eq. 15.21 a similar dipolar term is generated at the quartic term 
for rrii, a more detailed analysis is required to determine the fate of this quartic term. 

E. Isotropic triangular lattice 

This subsection will briefly comment on the modifications of our results for the case of 
the isotropic triangular lattice, with J' = J and full six-fold rotation symmetry. 

There is one more symmetry that needs to be considered: the 2n/3 rotation. Under this 
rotation, the visons of model II transform as 

R 27 r/3 : v x -> ^e-^V + -^=v 2 , v 2 -> —j=V! + -e^V (5.7) 

This PSG transformation is consistent with the enlarged U(l) gauge symmetry; while in 
model I, visons will be mixed with its complex conjugates, therefore on the isotropic trian- 
gular lattice only model II of visons is consistent. Further, the spinon minima are located at 
the commensurate wave vectors Q = (27r/3, 0) and — Q = — (27r/3,0), and therefore under 
translation spinons in model A will merely gain a phase factor, while in model B spinons 
will be mixed with their complex conjugates. Therefore on the isotropic triangular lattice, 
only model All is consistent with the enlarged U(l) gauge symmetries. 

In the mutual CS theory of model All on the isotropic lattice, g = 0, in the phase with 
both spinon and vison condensed, the VSB order parameter is described by the SU(2) vector 
v^a a v , which corresponds to degenerate stripe orders V s , Vy in Eq. (14.4D . and V z = v^o z v 
depicted in Fig. [3 these stripe orders are connected to each other through rotation R27r/3- 
Notice that, on the distorted triangular lattice, stripe order has the same symmetry as 
the square lattice VBS order. 

Because g = 0, the global symmetry of vison up to the forth order term is SU(2) in 
the mutual CS theory, and 0(4) in the theory with only visons. The GSM of the phase 
with both spinon and vison condensed is S 2 x S* 2 as far as the forth order terms are consid- 
ered. Therefore the transition between the Z 2 spin liquid and the VBS order is a 3D 0(4) 
transition^, and the transition between the nematic/VBS order to the spiral order in phase 
diagram Fig. [8] is a CP 1 transition. The PSG of visons allow a sixth order anisotropy term 
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on the isotropic triangular lattice^: 



L% = g&{v x vl + v 2 v\ + H.c.) 



(5.8) 



This term corresponds to the triple monopole process in the dual picture, which annihi- 
lates/creates three fluxes of gauge field a M . This triple monopole is expected to be relevant 
when gauge field a M is gapless, which will likely drive the transition between the nematic/VBS 
and the VBS phase to a first order transition. 

F. Multicritical point, s z = s v = 

We now study the multicritical point with both spinons and visons gapless: this is the 
point M in Fig. [31 

The most convenient way of studying M is likely via the non-CS theories in Eqs. (15. ip 
and (I5.2p . although these do not apply for the g > cases. This formulation should be 
amenable to direct numerical study. 

For analytic results, the only available tool is the 1/N expansion, for this we may as well 
work with the original CS theory in Eq. (II. ip . This expansion relies on the assumption that 
in Eq. (14.11) N ~ N v ~ k is large. The A term in Eq. 15.11 and 15.21 generalizes to terms with 
k powers of z a , and these are surely irrelevant for large k. So we ignore the influence of C x 
for models B in the 1/N expansion. 

A systematic 1/N expansion for the CP^ -1 -* model has been calculated previously^ 1 ^: 
the 1/N correction comes from the one-loop propagator of both the Lagrange multiplier A 
and gauge field a M : 



For instance, we can calculate the anomalous dimension t]n of gauge invariant operator 
z^T a z defined as A[z^T a z] = (D — 2 + rj N )/2, and T a is one of the generators of SU(N) 
algebra; note that this operator is the magnetic order parameter for model B, and the spin 
nematic order parameter for model A. In the CP^" 1 -* model, the anomalous dimension 77^ 



C = -\(d fl -ia li )z\ 2 + iX{\z\ 2 - !) + --■ 



(5.9) 



The one loop propagator of A and a M are: 




(5.10) 
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was calculated in detail in Reference |69_|, and the result is 

The second term on the r.h.s. of the equation above comes from the Lagrange multiplier, 
while the third term comes from the gauge field. 

After including the vison multiplet v a and mutual CS term, the A propagator is unaffected, 
while the gauge field propagator is modified: 

^J-(v-C^),n A ^( 1 + ^). (5J2) 



All the calculations can be carried out straightforwardly by replacing IT^ in Ref. |69| by 
ILa. When k ~ N, the correction from the vison and mutual CS theory to the anomalous 
dimension tjn is at the order of 1/N : 

32 128 1 
VN ~ 1 + 3^N ~ 3^N X 1 + 64F/(7r 2 iV 2 )' ^ 5 ' 13) 

All the other critical exponents can be calculated in a similar way. 



VI. BEYOND U(l) CHERN-SIMONS THEORIES 

In the previous sections, the global phase diagrams and nature of phase transitions were 
studied with the mutual CS theory. In this section we will try to look at these phases and 
transitions in a more intuitive and pictorial way. Let us start with the case with undistorted 
triangular lattice, which has a ground state with spiral antiferromagnetic order. This spiral 
order is described in terms of the z by Eq. (13.31) . 

Intuitively, to destroy magnetic order, the most straightforward way is to proliferate the 
topological defects in this magnetic order. The fundamental group of the GSM of spiral 
order is 7Ti[SO(3)] = Z 2 , which supports a topologically stable half vortex if rewritten in 
terms of z. These half vortices become full vortices of vectors n», i = 1,2,3 in Eq. ( 13. 31) . 
These vortices can be most easily understood as an ordinary vortices of two of the three 
vectors Ui, while keeping the third vector uniform in the whole 2d plane: see Fig. [TUJ The 
GSM of this state has isometry group SO (4). If the system has the enlarged 0(4) symmetry 
at the microscopic level, all the vortices with different uniform vector (UV) will have the 
same energy. However, the underlying symmetry is only SU(2) x PSG, so the energy of 
the vortices depends on the UV. The lattice symmetry guarantees that the vortices with 
different UVs have the same energy as long as the UVs can be transformed to each other 
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FIG. 10: Schematic of the orientation of the orthogonal vectors Hi, n 2 and n 3 around a vortex in 
the SO(3) GSM of the spiral antiferromagnet. One of the vectors has a constant orientation, while 
the other two precess by an angle of 2ir. 



through lattice symmetry transformations. There are in total three groups of vortices: 



All the flavors of vortices in each group have the same energy, while there is no symmetry 
to protect the degeneracy between different groups. For instance, UV n 3 can never be 
transformed in to n\ because of their opposite behavior under time reversal transformation. 
The first group of vortex has only one flavor, and if we only proliferate this vortex flavor, the 
spin orders of n\ and n 2 are destroyed, while the nematic order n 3 is preserved. This leads 
to a state with GSM S 2 , and it is the situation described by the mutual CS theory AI, All 
and AIII. The second and third groups have more than one flavor of vortices, with different 
flavors of vortices connected to each other through lattice symmetry transformations. If all 
the flavors of vortices in group 2 or 3 condense, the magnetic order is completely destroyed, 
and we expect these vortices to drive a direct transition from the v3 x antiferromagnetic 
spiral order to VBS order. The nature of this transition requires further study, and it cannot 
be naturally described by our mutual CS theory. In our theory, there are always two tuning 
parameters (s 2 and s v ), and it would require some fine-tuning to induce a direct transition. 
A theory based on a mutual Z 2 CS formalism has been proposed by another group^ to 
describe this direct transition. 

One can also condense one flavor of the second or the third group of vortices, which can 
be realized in the situation with strong repulsions between different flavors of vortices in one 
group. If the vortex with UV n*i is condensed, the vectors n 2 and n 3 are disordered, and the 



1, UV : n 3 



2, UV : ni 



U , \/3\ 
-2»i ± - 2-TI2, 



3, UV : n 2 



± — ni n 2 . 

2 2 



(6.1) 
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FIG. 11: The remnant spin order patten after proliferation of one single flavor of group 2 and 3 
vortices in Eq. (|6,ip . (a) Proliferation of the first flavor in group 2, the spin pattern is up-down- 
down, and the GSM is S 2 x Z3. Notice that the moment of the up spin site is twice as much as 
the down spin site, so the total magnetization is zero, (b) Proliferation of first flavor in group 3, 
the spin pattern is up-down-zero with GSM S 2 x Z3. 

remnant magnetic order is the up-down-down state in Fig. [TTb with zero total magnetization, 
and the GSM of this up-down-down state is S 2 x Z 3 . The S 2 corresponds to the direction 
of Hi, and the Z% corresponds to the choice of condensing the three flavors of vortices in 
the second group of Eq. (16. ip . which are connected to each other via translation along the x 
axis. If vortex with UV n 2 is condensed, the spin pattern becomes the up-down-zero state 
in Fig. [TTb . also with GSM S 2 x Z 3 . The transition driven by the condensation of vortices 
with UV rii and n 2 can no longer be described by the U(l) mutual CS theory, because in the 
phase with both spinons and visons condensed, the physical order parameter of the remnant 
spin order should be U(l) gauge invariant: thus now one needs a spinon z a such that the UV 
rii = z^a a z. However, under translation so-defined spinon z a becomes a linear combination 
between z a and e a/ gz^, which violates the U(l) gauge symmetry. 

To consistently describe the transition driven by proliferation of vortices with UV ni 
and n 2) a theory based on mutual Z 2 CS theory may be applicable, similar to the Z 2 CS 
formalism proposed for cuprates 3 - where the spinons and half- vortices of the superconducting 
phase are coupled together through mutual Z 2 CS fields. On the distorted triangular lattice 
antiferromagnets examined in this work, the mutual Z 2 gauge field also imposes the correct 
semionic statistics between the spinon and vison. However the Z 2 gauge field can only 
be conveniently formulated on the lattice, therefore we are unable to comment on all the 
universal properties of the transitions described by the mutual Z 2 CS theory. 

If the triangular lattice is distorted, the spin spiral state becomes incommensurate, and 
the vector n\ and n 2 can be transformed to each other via lattice translation. Therefore 
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there are only two groups of vortices: 



1, UV : n 3 , 

2, UV : Hi, rotation of Hi around H%. (6-2) 

The second group of vortices have an infinite number of flavors, and if one of these flavors 
proliferates, then the spin state becomes a collinear incommensurate spin density wave, 
which has GSM S 2 x S 1 . The S 2 corresponds to the remnant spin collinear order, while the 
S* 1 corresponds to translation along x of the incommensurate wave vector. 

VII. CONCLUSIONS AND EXPERIMENTAL IMPLICATIONS 

This paper has described examples of a general approach to describing the phases and 
quantum phase transitions of S = 1/2 antiferromagnets in two dimensions. Our examples 
were limited to models on the lattice of Fig. [TJ, because of its experimental importance. 
However, we expect that similar analysis should be useful e.g. on the kagome lattice. Our 
main results are summarized in phase diagrams, like that in Fig. [2j The same phase diagram 
was obtained earlier- in a more direct microscopic mean field theory, but the nature of the 
phase transitions and possible multicritical point was left open. Our present, 'dual' approach 
also immediately yields the required critical field theories. 

One of the transitions in our phase diagram in Fig. [2] is the Neel-VBS transition, which is 
described by a CP 1 field theory. This transition has been the focus of much recent numerical 
work 61 1 63 1 64 1 65 1 66 1 67. The numeric results on the S = 1/2 quantum antiferro magnet strongly 
support its effective description in terms of the CP 1 field theory. However, some results^ 1 ^ 
on system sizes larger than 50 x 50 indicate that the transition in the CP 1 model may well 
be weakly first-order. The multicritical point M has additional flavors of matter fields, and 
these make it less likely that M is first order. Consequently it would be useful to study M 
numerically using the theories in Eqs. ( 15. II) and ( 15.21) : this will help in describing M and 
the phase diagram in its immediate vicinity. This study can be done with or without the A 
coupling in Eqs. (15. ip and (15. 2p . 

As mentioned in Section [H a series of measurements on the distorted triangular lattice 
materials, X[Pd(dmit) 2 ]2, with different anisotropic coupling J'/ J reveal a possible direct 
transition between the Neel order and the VBS state, and one particular material with X = 
EtMesSb was suggested to be close to the quantum critical point. In our mutual CS theory, 
as well as the previously proposed deconfined critical point 3 -, this transition is described 
by the CP 1 model with irrelevant quadrupole operators. The quantum critical behaviors 
at finite temperature can be measured in material EtMe3Sb[Pd(dmit)2]2- For instance, the 
NMR relaxation rate 1/Ti scales as \jT\ ~ T w , where r/jv is the anomalous dimension of 
the Neel order parameter at the CP 1 fixed point. Even if the transition in the CP 1 model is 
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ultimately weakly first order, the critical regime could be observable at intermediate T, and 
compared to numerical estimates^ 3 - 1 ^ 1 ^ of t/jv- 

We conclude by commenting on the recent interpretation by Ref. |32| of experimental 
observations on the spin liquid compound k-(ET) 2 Cu2(CN) 3 . Ref. [32 suggested that the very 
low T nuclear magnetic resonance was controlled by the 0(4) criticality between the spiral 
and Z 2 spin liquid states, while the intermediate temperature nuclear magnetic resonance 
(NMR) could be modeled by a multicritical point where both the spinons and visons are 
gapless. Crucial to this interpretation was the requirement that the anomalous dimension 
of the magnetic order parameter, t]n, was smaller at the latter multicritical point than at 
the low T 0(4) critical point. 

Here, we have provided examples of such multicritical points, such as the point M in 
Fig. [2j For spinons in model B, the magnetic order parameter is z^T a z, and its anomalous 
dimension was computed in Section IV F\ where we found the result in Eq. (I5.13p . An 
important feature of this result is that the U(l) gauge fluctuations reduce the value of tjn 
from that obtained in the theory without the U(l) gauge field; the latter describes the 0(4) 
transition between the spiral antiferromagnet and the Z 2 spin liquid. Thus model B spinons 
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do fulfill the requirements for the experimental interpretation stated in Ref. 

On the other hand, in spinon model A, the order parameter z [ T a z represents the nematic 
order parameter, while the other gauge invariant spinon-monopole composite z t a v a a zM. h 
represents the spiral order parameter. In the case with large N v , the scaling dimension of 
Aib is expected to scale linearly with N v , and hence to systematically calculate the scaling 
dimension of z t a y <j a zA4b, we need an analytical technique beyond the 1/N expansion. 

In conclusion, we have proposed a unified theory in the mutual CS formalism to describe 
all the magnetic phases observed in a series of organic compounds, and discussed the phase 
transitions between these phases. Experimental implications for organic compounds in the 
X[Pd(dmit) 2 ]2 series and the k-(ET) 2 Z series were also discussed. 
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